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We investigate the ground state phase diagram of the S=l/2 two-leg XXZ spin lad- 
der system with an isotropic interchain coupling. In this model, there is the Berezinskii- 
Kosterlitz-Thouless transition which occurs at the XF-Haldane and the XF-rung singlet 
phase boundaries. It was difficult to determine the transition line using traditional meth- 
ods. We overcome this difficulty using the level spectroscopy method combined with the 
twisted boundary condition method, and we check the consistency. We find out that the 
phase boundary between XY phase and Haldane phase lies on the A = line. And we 
show that there exist two different XY phases, which we can distinguish investigating a XX 
correlation function. 

PACS numbers: 75.10.Jm, 75.40.Mg, 75.40. Cx 

I. INTRODUCTION 

Low dimensional quantum spin systems have attracted some attention. Especially, ground 
state properties are interesting since the quantum fluctuation often plays the dominant 
role. For example, S = 1/2 quantum spin ladder systems have been studied from both 
theoretical and experimental points of view 1 , related to Haldane's conjecture and high- 
T c superconductors. Haldane has predicted that the properties of integer spin chains are 
different from those of half-odd- integer spin chains 2 . Half odd integer spin chains do not have 
excitation energy gaps 3 . By contrast, for integer spin chains there is an energy gap between 
the ground state and the lowest excited state 2 . We can understand this fact, using the 
valence bond solid picture 4 . From the other point of view, Haldane's conjecture is related to 
spin ladder systems, since we can describe arbitrary spin-S" chains as S = 1/2 2S-\eg ladder 
systems approximately 5 . 

The model which we study in this paper consists of two coupled 5=1/2 XXZ chains. 
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We consider the following Hamiltonian, 
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(1) 



where a = 1, 2 are indices of chains (see Fig. 1), A is an anisotropy of the leg coupling, and 
J rung is the rung coupling. L is the system size (L = N/2, N is the number of sites). We 
shall discuss the boundary condition later. 

So far, S=l/2 two-leg quantum spin ladder systems have been studied by many au- 
thors. Strong and Mills have presented the phase diagram using the bosonization approach 6 . 
Watanabe et al. have discussed the same model using the bosonization with the aid of Wil- 
son's renormalization group method, and presented the phase diagram partially 7 . Narushima 
et al. have investigated Haldane-Neel transition 8 , using the density matrix renormalization 
group (DMRG) method 9 . These theoretical studies have been done mainly with antifer- 
romagnetic legs, since real materials have antiferromagnetic interaction, and those are in- 
vestigated in detail experimentally (e.g., Sri 4 Cu2404i). Sri 4 _ x Ca :r Cu2404i +< 5 is known as 
the superconductor under high pressure 10 . Recently, ferromagnetic legs cases attract much 
interest, for example, as Kolezhuk and Mikeska have discussed 11 . There exists a linear 
chain with ferromagnetic interaction, and we can expect the discovery of the ladder material 
with ferromagnetic interaction 12 ' 13 . Vekua et al. have studied ferromagnetic cases using 
the bosonization approach in the weak-coupling limit 14 . They found out that an extended 
gapless phase appears not only in the ferromagnetic rung coupling region but also in the 
antiferromagnetic rung coupling region. This phase was also referred by Strong and Mills 6 
too. The stability of this gapless phase for the A < region was discussed by Legeza and 
Solyom. 15 . A = 1 cases have been well studied using the bosonization 16,17 . In this case, the 
system is massive. Lecheminant and Orignac showed that the edge state does not exist in 
this model with J rung > under the open boundary condition (OBC) 17 . This means that 
there is no Haldane phase in J mng > region, since the edge state appears in Haldane phase 
under the OBC. 

There are several suggestions for the phase diagram of this model. However there is no 
unified view. There remain controversies in three points. The first point is the extent of the 
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XY phase. The second point is the possibility of existence of two different XY phases, and 
the third point is the discussion about the multicritical point near the Ferromagnetic phase. 

For the first point, there are some previous research. For example, Strong and Mills 6 
and Vekua et al. 14 have insisted that the XY phase extend to the A > region. Watanabe 
et al, who used the bosonization approach, have insisted the different extension in the 
A > region 7 . Narushima et al. 8 and Legeza et al. 15 have insisted that the XY phase does 
not extended in the A > region. They calculated the energy gap and the correlation 
function, using DMRG. But Narushima et al. and Legeza et al. did not determine the 
Berezinskii-Kosterlitz-Thouless (BKT) transition line, since the correlation length becomes 
large compared to the system size, especially near the XY phase. Also there are peculiar 
difficulties in determining the BKT transition point 18 . 

The XY phase belongs to the universality class of Tomonaga-Luttinger (TL) liquid 19 , 
which is characterized by a gapless excitation and a power-law decay of the correlation 
function. In general, one component TL liquid is described by the c = I conformal field 
theory (CFT) with the U(l) symmetry. One of the instabilities of the TL liquid is the 
BKT transition 20 . In the BKT transition, the traditional finite size scaling method (the 
phenomenological renormalization group 21 ) can not be applied 22,23 . Moreover, since there 
exists the logarithmic correction, it was difficult to determine the BKT transition point. 
Recently, combining the renormalization group with the symmetry consideration, the "level 
spectroscopy" method has been developed in order to overcome these difficulties 18 . Com- 
bining this method and twisted boundary condition 24 ' 25 , we will numerically determine the 
phase boundary between the XY phase and other massive phases in this paper. 

For the second point, the possibility of two different XY phases in ladder systems was 
pointed out by Vekua et al. 14 in a little different model. However, in their following paper, 
they did not comment on two different XY phases 26 . We discuss this point using a different 
approach in this paper. According to the Marshall-Lieb-Mattis's theorem, the XY\ phase 
is in J rung < region and the XY2 phase is in the J rung > region (see section. IV). 

The third point is the discussion about the multicritical point near the Ferromagnetic 
phase. There are different suggestions 6 ' 11 ' 14 . We clarified the slope toward the multicritical 
point at (A, J rung ) = (—1,0) concerning the contradicting results 11 ' 14 . 

We discuss the phase diagram of S=f/2 two-leg XX Z spin ladder systems in this pa- 
per (see Fig. 2) We organize this paper as follows. In the next section, we consider the 
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phase diagram under several limits. In section 3, we discuss the phase boundary, using 
the variational approach. In section 4, we discuss unitary transformations and correlation 
functions. In section 5, we determine some phase boundaries from exact results, considering 
the number of degeneracies of the ground state, and we consider the weak coupling region. 
In section 6, we numerically determine the phase boundary between the XY phase and 
the rung-singlet phase and between the XY phase and the Haldane phase, using the level 
spectroscopy method, and we discuss the multicriticality. The last section is the conclusion. 

II. SEVERAL LIMITS 

In this section, we consider the ground state in some limits. 

In the J rung = case, the system consists of two independent S — 1/2 XX Z chains, and 
it was solved exactly 27-29 . In the A < — 1 region, there is the ferromagnetic phase. In the 
— 1 < A < 1 region, the XY (spin fluid) phase appears. In the A > 1 region, the Neel phase 
appears. 

In the J rung — > — oo case, this model can be mapped onto an S — 1 XXZ spin chain. 
The S = 1 XXZ spin chain has been studied by many authors 30 ' 31 . In this system, there 
are four phases, i.e., the ferromagnetic phase, the XY phase, the Haldane phase, and the 
Neel phase. This Haldane phase is the gapped phase with the hidden Z 2 x Z 2 symmetry 
breaking. 

In the A — > — oo case, we can expect that ferromagnetic ordered phases appear in each 
leg directions. And we can classify them as follows. If J rU n g is positive, the ground state 
is called as the stripe ferromagnetic phase 14 . The staggered magnetization appears in the 
rung direction. This phase has the following correlation function: 

{SS^>« (-I)"*, (2) 

(see Fig. 3(a)). If J rung is negative, the ground state is the fully ferromagnetic phase. 

In the A — > +oo case, we can expect that antiferromagnetic (Neel) ordered phases ap- 
pear in each leg direction. If J rung is positive, the ground state is the Neel phase in the 
rung direction, (see Fig. 3(c)). If J rung is negative, the ground state is the ferromagnetic 
correlation in the rung direction. This phase is called as the stripe Neel phase. The stag- 
gered magnetization appears in the leg direction. This phase has the following correlation 
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function, 

(s^)«(-ir, (3) 

(see Fig. 3(b)). 

In the A ~ and J rung ^> 1 case, we can expect that the rung-singlet phase appears. 
This state is a nondegenerate, massive phase. We can understand this state as the direct 
product of singlets in the rung direction. 



III. VARIATIONAL APPROACH AND PHASE BOUNDARIES 

In this section, we roughly estimate phase boundaries using the variational method in 
some limits under the periodic boundary condition (PBC). 



A. Rung Singlet-Stripe Ferromagnetic phase transition 

At first, we consider the phase transition between the rung singlet phase and the stripe 
ferromagnetic phase. The variational energy of the pure rung-singlet state is given as follows: 

Ers _ (RS \H\ RS) ^3 
L ~ L - 4 J ™ng, 14 j 

where L is the system size (L = N/2, N is the number of sites). On the other hand, the 
variational energy of the pure stripe-ferromagnetic state is given as follows: 

E SF _ (SF|M|SF) _ 1 1 

~r _ i _ 2 4 Jmng - {b) 

So we can roughly estimate that the phase boundary between the rung singlet phase and 
the stripe ferromagnetic phase is A = — J rung in the A — > — oo and J rung — > oo limit. 

This phase transition is related to an spontaneously Z2 symmetry breaking. Therefore, 
this is the second order phase transition and its universality class may be the two dimensional 
(2D) Ising type. 



B. Neel-Rung Singlet phase transition 



Next, we consider the phase transition between the Neel phase and the rung singlet phase. 
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The variational energy of the pure Neel state is given as follows: 

E Neel _ (Neel Neel) 1 _ 1 
L ~ L ~ 2 4 rung ' 1 > 

So we can roughly estimate the phase boundary between the rung singlet phase and the Neel 

phase is A = J rung , in A — > oo and J rung — > oo limit. This phase transition is related to the 

Z2 symmetry breaking. Therefore, this phase transition is the second order phase transition 

and its universality class may be the 2D-Ising type. 

IV. UNITARY TRANSFORMATION AND PHASE TRANSITION 

In this section, we treat the bipartite systems. That is, the system size (L = N/2) is 
even under the periodic boundary condition, or the system size (L) is arbitrary under the 
open boundary condition, where N is the number of spins. 

A. XY phase - Ferromagnetic phase transition 

In the ferromagnetic coupling case (A < —1, J rung < 0), we found that the ferromagnetic 
phase appears. Vekua et al. have obtained that the phase boundary between XY phase and 
the ferromagnetic phase locates at A = — 1 using the instability analysis of the spin wave 
theory for the ferromagnetic phase 14 . Here we try to use another approach for this phase 
transition. We transform the original Hamiltonian (1) using a following unitary operator: 



U = exp 



(7) 



ITT 

(see Fig 4). This unitary operator transform S^ 2 j+i ~^ ~^a,2j+v Especially, in the A = — 1 
case, the transformed Hamiltonian is described as the pure Heisenberg model with a negative 
coupling constant. Now we can choose a set of basis vectors as eigenvectors which diagonalize 
S z . Then all off-diagonal elements of the transformed Hamiltonian have a negative sign. 
From Perron- Frobenius's theorem, the ground state is unique in the fixed = £V a J - 
space in the finite system, and we can choose that all the coefficients are positive. Since 
the transformed system has an SU(2) symmetry in the A = — 1 case, states with = 
0,±1,±2,--- are degenerate (SU(2) Ferro). This means that J rung < 0, A = — 1 line is 
the phase boundary, between the fully ferromagnetic phase (A < —1) and the XY phase 
(A>-1). 
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Especially, in the (J mng = 0, A = —1) case and the system size L(— N/2)), where L is 
even, there are the number of (L + l) 2 degenerated ground states, since this system has an 
577(2) x SU{2) symmetry. We can consider that this point is the multicritical point, where 
the BKT transition line meets the 2D Ising type phase transition line. 



B. Off critical case.l (J rU n g < 0) 

In this case, we consider the unitary transformation (7), (see Fig 4). After this unitary 
transformation, all off-diagonal elements of Hamiltonian become negative. From the discus- 
sion of the previous subsection, signs of correlation functions are represented in the original 
Hamiltonian (1) as follows, 

(• n' '(.v;:.,s; ; > >o. (8) 

This corresponds to the Ferromagnetic phase, the XY1 phase, Haldane phase and the Stripe 
Neel phase. 



C. Off critical case. 2 ( J rung > 0) 



In this case, we consider the following unitary transformation: 



U = exp 



ITi 



^ i S i,2j+l + S 2,2j) 



(9) 



S 2 2j ■ After 



(see Fig 5). This unitary operator transform S l 2 j+\ ~^ ~Si,2j+i an d S 2 2 j ~ 
this unitary transformation, all off-diagonal elements of Hamiltonian become negative. From 
the Marshall-Lieb-Mattis's theorem 32 , this case has following signs of the correlation function 
in the original Hamiltonian (1): 

(-1)^(-1)^(5^.>>0 (10) 

This corresponds to the Stripe Ferromagnetic phase, the XY2 phase, the rung singlet phase 
and Neel phase. 



D. Off critical case. 3 (A = -1, J rung > 0) 

In this subsection, we consider (A = — 1, J rung > 0) case. Using the unitary transfor- 
mation (7), we can see that the system has an SU(2) symmetry. From Mermin- Wagner's 
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theorem, since there is no long range order with the spontaneously continuous symmetry 
breaking, the ground state is not the Stripe Ferromagnetic phase. Since there is no soft 
mode with the wave number q = rr for the transformed Hamiltonian, there is no possibility 
that it is described as a k — 1,2 SU(2) Wess-Zumino-Witten model, which is massless. 

From the Marshall-Lieb-Mattis's theorem 32 , the sign of the correlation function trans- 
formed by unitary operator (7), 

{shst,) < 0. (11) 

Thus this ground state is not an SU (2) ferromagnetic phase. 

Therefore we can conclude that the ground state is the rung singlet phase in the (A = 
1, (/rung case. 

V. APPROACH FROM DECOUPLED CHAINS 

In this section, we discuss the phase diagram from the exact solutions. This system with 
/rung = consists of two independent chains. From the exact solution we know that S — 1/2 
XX Z chain with anisotropy |A| > 1 is massive, and the ground state is twofold degenerate. 
Now the system is two decoupled chains, therefore, the ground state is fourfold degenerate. 
In following subsections, we discuss the neighborhood of J rung = 0. 

A. Ferromagnetic phase - Stripe Ferromagnetic phase transition 

In A < — 1 region, we can find two different two ferromagnetic ordered phases, as sug- 
gested by Vekua et a/ 14 . The usual ferromagnetic ordered phase is two fold degenerate. 
However on the A < —1, J rung = line, the ground state is four fold degenerate. In the 
first case all spins of both chains are up; in the second case all spins of both chains are 
down; in the third case all spins of 1-chain is up and all spins of 2-chain is down; and in 
the forth case all spins of 1-chain is down and all spins of 2-chain is down (see Fig. (3) (a)). 
These four states have the same energy exactly on the J mng = 0. A nonzero 
breaks this degeneracy. This means that A < — 1, J run g = line is the phase boundary. In 
addition, the system, which consists of two independent gapped system is gapped. So this 
phase transition is the first order transition. 
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B. Neel phase - Stripe Neel phase transition 

The discussion in the previous subsection can be applied to the antiferromagnetic ordered 
phase similarly. 

Firstly, we consider the single S = 1/2 XXZ chain. From the exact solution, in A > 1 
region, the ground state is the antiferromagnetic (Neel) ordered state, which is two fold 
degenerate. Secondly, we consider independent double S — 1/2 XX Z chains. The ground 
state is fourfold degenerate in A > 1 region, since this total system is described as the 
direct product of the single chain system. These ground states are shown in Fig. 3 (b),(c). 
Lastly, we consider the coupled S = 1/2 XXZ chains described by the Hamiltonian (1). 
The rung-coupling term breaks the four fold degeneracy in the above case. In J rung < 
case, the ferromagnetic ordered state in the rung direction (see, Fig. 3 (b)), which is called 
as the stripe Neel phase. In J rung > case, the antiferromagnetic ordered state in the rung 
direction (see, Fig. 3 (c)), which is called as the Neel phase. Both states are antiferromagnetic 
ordered in the leg direction. 

Since the system is gapped in the J rung = case, this phase transition between the Neel 
phase and the stripe Neel phase is the first order transition. 

C. Weak coupling region 

In this subsection, we consider the weak coupling S=l/2 XXZ chain, — 1 < A < 1 and 

| J rung | 1 • 

At first, we consider two independent chains, J rU n g = case. This case is exactly solved 
using Bethe ansatz 27-29 . One chain in the — 1 < A < 1 region is described by TL liquid 19 . 
Two independent chains are described by the direct product of two TL liquids. This model 
has critical properties in extended regions in the parameter space. Here we consider the case 
that there are some perturbations for this system. 

These systems can be analyzed using the bosonization method 6 ' 14 . Then two decoupled 
chains are described in two bosonic fields, (0i,#i) and (02, #2) of each other. We obey the 
notation used in Ref. 14 and introduce the symmetric and antisymmetric combination of the 
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bosonic field, 



<f>± OC 01 ± 02, 

e± oc e 1 ± 2 . 



(12) 



We consider the case with J rung ^ 0. Behavior of the symmetric field is governed by the 



effective sine-Goldon model, cos y i s always relevant, cos y / 2nK + (f) + is relevant for 

A > 0, irrelevant for A < 0, where K± are constant 14 . 

The phase transition between the XY phase and the nondegenerate massive phase is the 
K = 4 BKT type transition when the system has a simple U(l) symmetry, where K is the 
TL parameter. In this case, a = 0,q = ir) mode is always massive, and a (S^ — 0, q — 0) 
mode becomes massless. 

We can distinguish two different XY phases. One XY phase (XY1) is in J rung < 0, 
another XY phase (XY2) is in J rung > region. In the J rung > region, from the 
Marshall-Lieb-Mattis's theorem, (Sf^S^j) < 0. On the other hand, in the J mng < re- 
gion, (S*jS%j) > 0. Therefore these two XY phases have different symmetry. A J rung = 
line is the second order phase transition line. The system on the J rU n g = line is described 
as two component TL liquid. This second order phase transition is explained as that one 
component of the TL liquid remains massless, but another component become massive. 

Here we summarize this subsection in the CFT language. In the |A| < 1 region, the 
system is described by the central charge c = 2 CFT on the J run g = line. And there exists 
a region which is described as c = 1 CFT in the < | J rung | 1 and — 1 < A < 0. 



In this section, we determine the critical points and the universality class. We show the 
numerical results of iV = 12, 16, 20, 24 systems using the exact diagonalization, and the level 
spectroscopy method with the twisted boundary condition 18 ' 24 ' 25 . 




VI. PHASE TRANSITIONS AND NUMERICAL RESULTS 



A. XY-Rung singlet transition 



We can expect that this phase transition is the BKT type. Using the level spectroscopy 
method, we determine critical points of this phase transition. Considering symmetries, we 
can identify critical points as cross points of two low-lying excitations. One has quantum 
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numbers Sj. = ±2, P = 1, g = under the periodic boundary condition (PBC), another has 
quantum numbers — 0, P* = 1 under the twisted boundary condition (TBC) (5^ = 
—S*f +L ). Here is total magnetizations of the system = ^^y, and P is the 
parity defined under the following transformation (Sfj <-> S'fi-j+i), P* is the parity under 
the TBC. g is the wave number defined under the periodic boundary condition, L is the 
system size which is a half of the number of sites. We show crossing points in Fig. 6. 

However values of crossing points have size dependence. We can remove logarithmic 
corrections, which is proportional to 1/ InL, using the level spectroscopy method. However, 
there remain other corrections which come from the lattice structure. For example, there 
are correction terms from the x = 4 irrelevant field (L_2-L_2l) 33 . Therefore we extrapolate 
crossing points as follows, 

•/rung (L) = (oo) + a± + o (j^ (13) 



or 

A cross (L) = A cross ^ + + O (j^j (14) 

where a is a fitting parameter. We neglect higher order terms o (^) . In Fig. 7, we show the 
size dependence of the crossing points and its extrapolation. Then we can find the phase 
boundary under the infinite system limit in Fig. (8). 

In order to confirm the consistency of our results, we need check the universality class 
using CFT 34 . In this case, the critical theory which describes the transition line is the c = 1 
CFT. The central charge appears as an universal finite size correction for the ground state 
energy under the PBC 35-37 , 

irv 

where L is the system size, e is the ground state energy per site in the infinite system, v is 
the spin wave velocity. Now o ^ ^^3 ) correction terms are small enough for numerical data, 
so we neglect it. In order to numerically determine the central charge, we should obtain v 
in addition to the ground state energy. We can obtain the spin wave velocity as follows: 

«< £ > = sM« = t)- b ')- (16) 

where q is a wave number. Then we extrapolate v (L) as 

v (L) — v (oo) + a-^r + b— + higher order, (17) 



E g (L)=eL--(c + o(— 3 ) ) (15) 
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TABLE I: P = 1* is even parity under TBC, P = -1* is odd parity under TBC, , q is the wave number, x 
is the scaling dimension, yi oc 1/lnL 
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2 4 
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2^4 2 


T TBC 
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TBC 


1 1 Vl I V2 

2 ~ 4 ~ 2 


T TBC 


±4 


1 





PBC 
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£±4,0 





1 





PBC 


2-yi(i + f) 


%marg 





-1 





PBC 


2 + yi 







1 





PBC 


2 + 2 yi (l + f) 





where a and b are fitting parameters. We neglect higher order terms. 

Here we define the effective central charge in numerical calculations as follows, 

E g (L) = eL-^, (18) 

where c is the effective central charge. The effective central charge is equivalent to the central 
charge on the critical point or in the TL liquid phase. And this changes rapidly from c = 1 
(the TL phase) to c — > (the massive phase) 23 . In Fig. (9), we show the effective central 
charge obtained on the transition points. In this figure, we find that the central charge 
decreases from c = 2 (J rung = 0) to c = 1 (far from J rung = 0). This reflects Zamolodchikov's 
c-theorem 38 . And we consider the phase transition between the rung-singlet phase (massive 
)and the XY phase (massless). We see that the effective central charge c is 1 in normal XY 
phases, and c is 2 in XY phase in two independent chains J rung = 0, and c becomes zero 
increasing the system size (L — > oo) in the rung-singlet phase (see Fig. (10)). 

And then we should calculate scaling dimensions removing logarithmic corrections. From 
CFT, scaling dimensions in the finite size system under the PBC are related to excitation 
energies as follows 39 : 

Ei (L)-E g (L) = *?p, (19) 

where Ei is an excitation energy, v is the spin wave velocity, L is the system size, Xi is the 
scaling dimension. In fact, there exists additional logarithmic corrections. From table (I), 
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we obtain the following relation: 

2x ±2fi + x^+x^ s = 2, (20) 

removing leading logarithmic corrections. In Fig. 11, we show its size dependence. 

Further more, we check the consistency of BKT phase transition. We calculate the ratio 
of two excitation energies, — 1 and = 2. 

AE 1 = E(S Z T = 1) - E g , 
AE 2 = E(S Z T = 2) - E g , 

' - zf • < 21 > 

where E s is the ground state energy, E(Sj, = 1, 2) is the lowest energy with the quantum 
number = 1, 2. This / is equivalent to the ratio of scaling dimensions on the TL liquid. 
For Gaussian model 40 , scaling dimensions are described as 

1 fm 2 \ 
x m , n = -l — +n 2 Kj, (22) 

under the PBC, where m is related to the magnetization of the system, n is not related to 
any conserved quantities, K is the TL parameter, m and n are integers. From eq.(21) and 
eq. (22) / should be 4 in one component TL liquid phase. 

However, when the system consists of two independent chain, this / should be 2 in TL 
liquid phase for the following reason. In S z = 1 case, one chain has the ground state energy 
E g , another chain has the excitation energy E(S Z = 1), thus the system has the excitation 
energy AEi. In S z = 2 case, each chain has the excitation energy E(S Z = 1), so the system 
has the excitation energy 2A£ , i (see Appendix A). In this case, / should be 2. 

On the other hand, in a massive phase, the energy of two independent magnons is twice 
the energy of one magnon in L — > oo limit. If the interaction of magnons is repulsive, then 
the ratio is more than 2 in the finite system. If the interaction of magnons is attractive 
then the ratio is less than 2 in the finite system. We show the ratio / in Fig. 12. This 
figure support the second order phase transition between the XY1 phase and XY2 phase, 
together with Fig. 10. 
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B. XY -Haldane transition 

The procedure in the previous subsection can be applied to the XF-Haldane transition, 
since we can expect that XF-Haldane transition is BKT type. In this case, we can identify 
critical points as cross points of two low-lying excitations. One has quantum number = 
±2, P — 1, q — under the PBC, another has quantum number = 0, P* = —1 under the 
TBC. We show crossing points in Fig. (13). We can find the transition line is on A = 0. As 
before, we calculate the effective central charge and the scaling dimensions after removing 
logarithmic correction. We show results in Fig. 14 for the central charge, and in Fig. 15 for 
the scaling dimensions. As a result, we obtain the phase boundary is on a A = line. This 
is analogous to the phase boundary between XY phase and Haldane phase in S — 1 XX Z 
spin chain 30 ' 31 . 

We have found analytically that S=l XY spin chain has an additional SU(2) symmetry 
under the OBC and an artificial boundary condition(ABC) in our previous paper 41 . In this 
paper, we show that there is an additional SU(2) symmetry for S = 1/2 two-leg XX Z 
spin ladder system with A = under the OBC and ABC (see Appendix B and Ref. 41 ). 
This supports our numerical calculations, since an SU (2) symmetry is related to the BKT 
transition 42 . 

C. Multicritical point (A, J rung ) = (0,0) 

In this subsection, we investigate excitations near the multi-critical point J rung = and 
A = 0. On this point, BKT transition line meets Gaussian transition line. BKT transition 
line is a phase boundary between XY phase and Haldane phase, and between XY phase 
and rung-singlet phase. Gaussian transition line is a phase boundary between rung-singlet 
phase and Haldane phase. BKT transition line consists of two parts, one is the crossing line 
(S% = ±2,P = l,q = under the PBC) and (S%. = 0, P* = 1 under the TBC), another is 
the crossing line (Sj< = ±2, P = 1, q = under the PBC) and (S^ = 0, P* = —1 under the 
TBC). 

We show excitation energies with the following quantum numbers and boundary condi- 
tions: 



• S% = ±2,P = l,q = Q under the PBC 
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• = 0, P* = 1 under the TBC 

• = 0, P* = -1 under the TBC 

in Fig. (16) along the J run g = A line and in Fig. (17) along the J rU n g = — 4A line in the 
parameter space. We show these crossing points in Fig. (18). This show that (A, J mng ) = 
(0, 0) is a multicritical point. Wl-Haldane transition line and W2-rung singlet transition 
line are continuously connected but these are not smoothly connected. 

From bosonization studies 6 ' 14 , W-rung singlet transition line seems to smoothly connect 
W-Haldane transition line. However, this is not consistent with CFT analysis. Our nu- 
merical calculations support CFT analysis. We can consider that this reflects an additional 
SU(2) symmetry on A = line 41 . 

D. Multicritical point (A, J rung ) = (-1,0) 

In this subsection, we discuss the other multicritical point (A, J rung ) = (—1, 0). This mul- 
ticritical point is among the XY phase, the Ferromagnetic phase, the Stripe Ferromagnetic 
phase, and the rung-single phase. This multicritical point has been studied by Kolezhuk 
and Mikeska 11 and Vekua et al. u . Kolezhuk and Mikeska discussed this, mapping onto 
the anisotropic non-linear a model. The renormalization group analysis of the anisotropic 
non-linear a model was developed by Nelson and Pelcovits 43 . Vekua et al. also discussed 
the ferromagnetic leg case, using the bosonization 14 . Outlines of the phase diagram are the 
same, however there is different point. Kolezhuk and Mikeska suggested that the shape of 
the phase boundary between XY phase and rung-singlet phase is the exponential function 
Jrung oc exp (—a/ | A + 1|), where a is a positive constant. On the other hand, Vekua et al. 
suggested that the shape of phase boundary between the XY phase and rung-singlet phase 
is the linear line on A — J rung plane. Our numerical result strongly supports Kolezhuk's sug- 
gestion (See Fig.8). This shows that non-linear a model with anisotropy is more appropriate 
to describe the region near the multicritical point. 

VII. CONCLUSIONS 

In this paper, we have studied the ground state phase diagram of S=l/2 two-leg spin 
ladder system. We have accurately numerically determined the phase boundary between 
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the XY phase and the Haldane phase, and between the XY phase and the rung-singlet 
phase, analyzing the exact diagonalization data using the level spectroscopy method, TBC, 
and CFT. And we have checked the universality class. As a result, in — 1 < A < region, 
the XY phase extends over the J rung = line from the J rung < region to the J rung > 
region. And there does not exist the XY phase in A > region except on the J rung = line, 
since the BKT transition line is the A = line. We can understand this result, considering 
an additional SU(2) symmetry 41 . 

We have roughly discussed phase boundaries between the rung-singlet phase and Neel 
phase, and between the rung-singlet phase and the stripe ferromagnetic phase, using the 
variational method. We think that these phase transitions are the second order phase 
transitions, considering the symmetry, and that the universality class is the Ising type. 

We have determined phase boundaries between the stripe ferromagnetic phase and the 
fully ferromagnetic phase, and between the Neel phase and the stripe Neel phase, considering 
degeneracies of ground states. These phase transitions are of the first order type. 

We find that there are two different XY phase. We can distinguish these XY phases, 
considering XX (or YY)correlation function. This is based on the Marshall-Lieb-Mattis's 
theorem. One XY phase (XY2) is the J mng > region, other XY phase (XY1) is the J run g < 
region, and it is understood that the second order phase transition occurs numerically 
between the XY1 phase and the XY2 phase from Fig. 10 and Fig. 12. 

J rung = 0, A = — 1 point is a multicritical point, among the rung-singlet phase, the stripe 
ferromagnetic phase, the fully ferromagnetic phase, and the XY phase. Seeing the J rung > 
region, the phase diagram is quite characteristic. This phase diagram is similar to one of 
the anisotropic nonlinear a model without the topological term whose phase diagram and 
the renormalization group flow have been discussed by Nelson and Pelcovits 43 . The same 
analysis can be applied to this problem. 
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APPENDIX A: EXCITATIONS OF THE TWO INDEPENDENT CHAINS SYSTEM 

In S — 1/2 two-leg ladder system, the J rung = case is described as two independent 
S = 1/2 XXZ chains. An S — 1/2 XXZ chain is exactly solved 28 , using the Bethe ansatz 27 . 
This model is described by the Hamiltonian 

H = JJ2 (^m + S * S li + AS ? S Ui) ■ (Al) 

i 

For A > 1, the system becomes the Neel (Anti Ferromagnetic) phase, whereas for A < — 1, 
the system becomes the Ferromagnetic phase. In the intermediate case — 1 < A < 1, the 
system becomes XY (spin fluid) phase. 

In this appendix, we discuss the excitation structure of two independent S=l/2 XY spin 
chain. 

In XY phase, the system is described as c = 1 gaussian model. Now we have two 
independent chains. We describe the state of all systems as We think that is the 
state of chain 1, and \1> 2 is that of chain 2. 

• vj> = \& 1 i]/ 2 

• Q — Qi + <?2 : wave number 

• Sj — Sf 1 + Sf 2 : magnetization 
Excitations of the single chain are, for example, 



• 


q 


= o,sJ 


= 0,x = 


: the ground state 


• 


Q 




= 0,x = 


1 


• 


q 


= o,s? 


= 0,x = 


2 


• 


q 


= ir,Sj 


= ±l,x 


l 

4' 



Excitation of the system are as follows, 



• (q = ir, Sj 1 = 0, x = 1) ® (q = ir, Sj 2 = 0, x = l) = (q = 0(2tt), Sj = 0, x = 2) 

• (q = 0, S^ 1 = 0, x = 2) (g) (q = 0, S? 2 = 0, x = 0) = (q = 0(2tt), Sj = 0, x = 2) 

. (q = tt, SP = +l,x = \)®(q = tt, S T Z 2 = -\,x=\) = (q = 0(2vr), = 0, x = \) 



18 



The last one transforms as follows, under the spin reversion S z <-> —S z , 
(q = 7T, SJ = +l,x = \)®(q = tt, S T Z =-l,x = \) 

(q = 7T, S Z T = = i) <g> (g = 7T, Sj = +l,x = ±) 

This has the odd rung-parity. This excitation is relevant, it does not appear in the space 
where the rung-parity is restricted to even. 

APPENDIX B: AN ADDITIONAL SU(2) SYMMETRY OF THE S=l/2 TWO-LEG SPIN 

LADDER SYSTEM 41 

We present a brief review of an additional SU(2) symmetry of the S — 1/2 two-leg ladder 
system. In the 5 = 1 XY chain and the S — 1/2 two-leg ladder system, there exists an 
additional SU (2) symmetry besides the usual spin SU (2) symmetry. In our previous paper 41 , 
we show that the one dimensional spin-1 XY model has an additional SU (2) symmetry. And 
we also show the case of S — 1/2 two- leg ladder system. 

We consider the Hamiltonian (1). At first, we introduce the following local operators: 

*j=StA ^ = \(^ J + S^). (Bl) 

After the simple calculation, we obtain the commutation relation of SU(2), 

[s*, Sf] = ±5 jk sf, [S+ 5~] = 25 jk s* (B2) 

The operator £\ §j commutes with the Hamiltonian (1), but operators £\ sf do not. 
Thus we introduce new non-local operators, 

4 = S tj S t> I[ H S h S h) *j = \ + S Q ■ (B3) 
1=1 

New operators satisfy the commutation relation of SU(2), 

[s*, 4] = ±5 jk sf [ a + sj;] = 25 jk s* (B4) 

and 
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where = and Sj are total operators. 

After the some calculation, we can show that these total operators commute with the 
Hamiltonian (1) with A = under the open boundary condition and an artificial boundary 
condition. For a discussion in detail, please see our previous paper 41 . 
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FIG. 1: Schematic structure of a S=l/2 two-leg spin ladder of Eq (1) 
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FIG. 2: Schematic phase diagram in the S=l/2 XXZ ladder system on A — J rU n g plane. Dotted lines are 
phase boundaries. 
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FIG. 3: (a) Schematic picture of the Stripe Ferromagnetic state. This phase has ferromagnetic order in 
the leg direction, and antiferromagnetic order in the rung direction, (b) Schematic picture of Stripe Neel 
state. This phase has antiferromagnetic order in the leg direction, and ferromagnetic order in the rung 
direction, (c) Schematic picture of Neel state. This phase has antiferromagnetic order in the leg direction, 
and antiferromagnetic order in the rung direction. 



Sl,n-1 Sl,n Sl,n+1 Sl,n+2 Sl,n+3 

— o — • — o— 



— o — • — o— 

S2,n-1 S2,n S2,n+1 S2,n+2 S2,n+3 



FIG. 4: Schematic figure for the unitary transformation in the J rU n g < case. Empty circles are uncchanged. 
Full circles are S± d -» -S^ 
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FIG. 5: Schematic figure for the unitary transformation in the J rU n g > case. Empty circles are uncchanged 



Full circles are Sgj -> -S*,- 
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FIG. 6: Excitation energies of L = 12(= N/2), A = —0.5 near the XY-Rung Singlet transition point, x's 
are iSf. = ±2,g = 0,P = 1 under the periodic boundary condition. +'s are — 0,P* — 1 under the twisted 
boundary condition. *'s are Sj, = 0,P* = —1 under the twisted boundary condition. 
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FIG. 7: Size dependence of the crossing points. The extrapolated value is 0.3727 at A = —0.5. Dotted line 
is determined using a linear least method. 
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FIG. 9: The effective central charge on XY-Rung Singlet transition line, x's are extrapolated using L(= 
N/2) = 6, 8, 10. +'s are extrapolated using L(= N/2) = 8, 10, 12. 




FIG. 10: The effective central charge c as a function of J run g at A = —0.5. x are extrapolated using 
N = 12, 16, 20. + are extrapolated using N = 16, 20, 24. 
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FIG. 11: Size dependence of scaling dimensions removed logarithmic correction on the XY-Rung Singlet 
transition line for the system size N = (2L) = 12, 16, 20, 24. 
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FIG. 12: The ratio of energy gaps as a function of J run g at A = —0.5 and its size dependence. 
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FIG. 13: Excitation energies of L = 12(= N/2), J rung = —0.5 near the XY-Haldane transition point, x's 
are Sf. = ±2,g = 0,P = 1 under the periodic boundary condition. +'s are Sj. — 0,P* — 1 under the twisted 
boundary condition. *'s are 5f. = 0,P* = —1 under the twisted boundary condition. 
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FIG. 14: The effective central charge on XY-Haldane transition line, x 's are extrapolated using L(= N/2) = 
6, 8, 10. +'s are extrapolated using L(= N/2) = 8, 10, 12. 
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FIG. 15: Size dependence of scaling dimensions removed logarithmic correction on the XY-Haldane transition 
line (A = 0) for the system size N(= 2L) = 12, 16, 20, 24. 




FIG. 16: Excitation energies of L = 12(= N/2), on A = J rung line near the multicritical point, x's are 
5f, = ±2,q = 0,P = 1 under the periodic boundary condition. +'s are Sf, = 0,P* = 1 under the twisted 
boundary condition. *'s are = 0,P* = —1 under the twisted boundary condition. 
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FIG. 17: Excitation energies of L = 12(= N/2), on A = — 4J rung line near the multicritical point, x's are 
Sj, — ±2,q = 0,P = 1 under the periodic boundary condition. +'s are Sf. = 0,P* = 1 under the twisted 
boundary condition. *'s are Sj. = 0,P* = —1 under the twisted boundary condition. 
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FIG. 18: Extrapolated cross points of excitation energies, x's are cross points of the excitation with 
Si? = ±2,g = 0,-P = 1 under the periodic boundary condition and the excitation with = 0,P* = 1 under 
the twisted boundary condition. +'s are cross points of the excitation with Sj- = ±2,g = 0,P = 1 under the 
periodic boundary condition and the excitation = 0,P* = — 1 under the twisted boundary condition. 



